Goal-Oriented Adaptivity and Multilevel Preconditioning for the Poisson-Boltzmann Equation
نویسندگان
چکیده
In this article, we develop goal-oriented error indicators to drive adaptive refinement algorithms for the Poisson-Boltzmann equation. Empirical results for the solvation free energy linear functional demonstrate that goal-oriented indicators are not sufficient on their own to lead to a superior refinement algorithm. To remedy this, we propose a problem-specific marking strategy using the solvation free energy computed from the solution of the linear regularized Poisson-Boltzmann equation. The convergence of the solvation free energy using this marking strategy, combined with goal-oriented refinement, compares favorably to adaptive methods using an energy-based error indicator. Due to the use of adaptive mesh refinement, it is critical to use multilevel preconditioning in order to maintain optimal computational complexity. We use variants of the classical multigrid method, which can be viewed as generalizations of the hierarchical basis multigrid and Bramble-Pasciak-Xu (BPX) preconditioners.
منابع مشابه
Electrostatic analysis of the charged surface in a solution via the finite element method: The Poisson-Boltzmann theory
Electrostatic potential as well as the local volume charge density are computed for a macromolecule by solving the Poisson-Boltzmann equation (PBE) using the finite element method (FEM). As a verification, our numerical results for a one dimensional PBE, which corresponds to an infinite-length macromolecule, are compared with the existing analytical solution and good agreement is found. As a ma...
متن کاملThe adaptive multilevel finite element solution of the Poisson-Boltzmann equation on massively parallel computers
Using new methods for the parallel solution of elliptic partial differential equations, the teraflops computing power of massively parallel computers can be leveraged to perform electrostatic calculations on large biological systems. This paper describes the adaptive multilevel finite element solution of the Poisson-Boltzmann equation for a microtubule on the NPACI IBM Blue Horizon supercompute...
متن کاملAdaptive multilevel finite element solution of the Poisson-Boltzmann equation II. Refinement at solvent-accessible surfaces in biomolecular systems
We apply the adaptive multilevel finite element techniques described in [20] to the nonlinear Poisson-Boltzmann equation (PBE) in the context of biomolecules. Fast and accurate numerical solution of the PBE in this setting is usually difficult to accomplish due to presence of discontinuous coefficients, delta functions, three spatial dimensions, unbounded domains, and rapid (exponential) nonlin...
متن کاملThe Adaptive Multilevel Finite Element Dolution of the Poisson-Boltzmann Equation on Massively Parallel Computers
By using new methods for the parallel solution of elliptic partial differential equations, the teraflops computing power of massively parallel computers can be leveraged to perform electrostatic calculations on large biological systems. This paper describes the adaptive multilevel finite element solution of the Poisson–Boltzmann equation for a microtubule on the NPACI Blue Horizon— a massively ...
متن کاملExtensions to Study Electrochemical Interfaces - A Contribution to the Theory of Ions
In the present study an alternative model allows the extension of the Debye-Hückel Theory (DHT) considering time dependence explicitly. From the Electro-Quasistatic approach (EQS) done in earlier studies time dependent potentials are suitable to describe several phenomena especially conducting media as well as the behaviour of charged particles in arbitrary solutions acting as electrolytes. Thi...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- J. Sci. Comput.
دوره 52 شماره
صفحات -
تاریخ انتشار 2012